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Charge transport in GaN quantum well (QW) devices grown in the non-polar direction is theoretically investi- 
gated . Emergence of a new form of anisotropic line charge scattering mechanism originating from anisotropic 
rough surface morphology in conjunction with in-plane built-in polarization is proposed. It is shown that 
such scattering leads to a large anisotropy in carrier transport properties, which is partially reduced by strong 
isotropic optical phonon scattering. 



7 The growth of non-polar m-plane (llOO) and a-plane 
s (1120) GaN has attracted a lot of interest recently 1 ' 2 . 

9 Though the built-in polarization field in traditional 

10 polar-GaN (c-plane) has been exploited to achieve dopant 

11 free high electron mobility transistors (HEMTs), for op- 

12 tical devices, polarization field plays a negative role due 

13 to quantum confined Stark effect (QCSE). Moreover, c- 

14 plane GaN-based enhancement-mode (E-mode) HEMTs 
is have very low threshold voltage (V t h ~ 1 V) due to 

16 inherent presence of polarization induced electron gas 

17 (2DEG). High threshold voltage (Vth ~ 3 V) is required 
is for high- voltage switching operations. So the recent trend 

19 is to explore optical and transport properties of GaN 

20 grown in the non-polar direction. Recently E-mode tran- 

21 sistor with high threshold voltage (> 2 V) has been 

22 achieved 3 with non-polar GaN hetero junction. 

23 A characteristic feature of non-polar GaN surfaces grown 

24 directly on foreign substrates is extended stripe or slate 

25 like morphology perpendicular to the c axis 2 ' 7 10 , regard- 

26 less of growth methods. The origin of this rough sur- 

27 face morphology has been attributed to i) replication of 

28 substrate morphology 4 , ii) extended basal-plane stack- 

29 ing faults (BSFs) 2 ' 7 , and hi) anisotropic diffusion bar- 
so rier of Ga adatoms 5 ' 6 . This striated morphology has 

31 been conjectured qualitatively to be responsible for ex- 

32 perimentally observed conductivity anisotropy for bulk 

33 GaN films 11 as well as thin GaN-QWs 12 . No microscopic 

34 theory is available in existing literature for a quantitative 

35 estimation of this electrical anisotropy. This paper devel- 

36 ops a microscopic theory of anisotropic carrier transport 

37 in non-polar GaN QWs with a quantitative estimation of 

38 the transport anisotropy. 

so Let us consider a thin non-polar GaN QW of thickness 

41 a sandwiched between two aluminum nitride (A1N) lay- 

42 ers as shown in Fig. la). A common source of disorder 

43 for non-polar GaN QW is striped surface morphology 

44 as mentioned in the previous section. We model this 

45 rough surface morphology as a local variation of QW 

46 thickness ( see Fig.l) The variation of QW thickness 

47 alone causes local shifts of the conduction band edge, 
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FIG. 1. Non-polar GaN quantum well: a)QW of width a 
sandwiched between A1N barriers, b) schematic diagram of 
interface roughness at GaN/AIN interface, c) polarization in- 
duced line charge at each step edge of roughness. 

48 resulting in carrier scattering commonly known as in- 

49 terface roughness scattering (IR) in literature. The in- 

50 terface roughness (IR) can be modeled by local thick- 
si ness fluctuations A(x) of the non-polar GaN QW with 

52 a spatial correlation (A(x)A(O)) = A 2 exp[— |x|/\/2A], 

53 where A is the average height of roughness and A is 

54 the in-plane correlation length between two roughness 

55 steps. Denoting the envelope function of conduction elec- 

56 trons in the n th subband of GaN QWs as |n,r, k) = 

57 y/2/a siri(7rnz / a /)e lk ' r / y/L x L y , the square of the un- 
ss screened intra-subband IR matrix element of scattering 

59 from initial momentum state ki(fc$,0) to a final momen- 

60 turn state kf(fc/,0') in the m th subband can be written 
ei as 

2 m 2 ( Tr 2 h 2 A ^ 2V2A 
\v IR (q x )\ = —{-^) 2+(&A)2 ^.o. 

62 where, q = kf — ki ,L X and L y are the macroscopic 

63 lengths of the QW in x and y directions, is Kronecker 

64 delta function and Ti is the reduced Planck constant. The 
es problem of IR scattering is well known and has been in- 
ee vestigated for two and one dimensional electron gases in 
67 many semiconductors, including nitrides. What is new 
es in non-polar structure is the polarization induced bound 
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100 the matrix element of the dipole potential summed over 

101 all scatterers is given by 



FIG. 2. schematic diagram of momentum change in 
anisotropic elastic scattering. a):momentum space. Clearly 
transverse wave vector (k y ) of carriers parallel to line charges 
is conserved while the net momentum change is |2fccos0|, b) 
in real space electron's movement is shown. 



charges associated with each interface roughness step. 
The thickness modulation of QW leads to GaN/AIN het- 
erojunctions at each roughness center as depicted in Fig.l 
b). The difference of in-plane polarization of GaN and 
A1N creates line-charge dipoles by inducing fixed charges 
at opposite faces of each step (as shown in Fig.l.b). As- 
suming each rough step is infinitely extended along the y 
direction (see ref. 6-9), the electrostatic potential at any 
point r (x, z) in the QW arising from a polarization 
line-charge dipole of a roughness step at a point (a^,0) 
of height A(xi) and lateral width d is given by 



V(x,y,z) = v+(x,y,z) - v-(x,y,z) 



In 
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so where e is the electron charge, eo is the free- space per- 
8i mittivity, n is the relative dielectric constant of GaN 



82 and X n (xi) 



A(xi)/e is the effective 
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83 line charge density. Note that, the potential is in- 

84 dependent of y due the symmetry of the problem, 
ss The scattering matrix element of transition from state 
se \n,r,ki) to state \m,r,kf) can be written as v(xi,q) — 
87 (m, r, kf\v(x, y, z)\n, r, where 



v(xi,q) = ( 
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88 In the above equation, F nm (q x a) is the form factor aris- 

89 ing from the quasi-2D nature of the electron gas. For our 

90 choice of the envelope function, the form factor can be 

91 calculated analytically. It approaches unity (F nm (q x a) —> 

92 1), both for long wavelength (q — ^ 0) as well as for for a 

93 very thin QW (a — )> 0). 

9i Eq. (3) represents the Fourier-transformed electrostatic 

96 potential from a line charge dipole associated with a sin- 

97 gle roughness step. As two steps are correlated, the 

98 dipole potential arising from them are also correlated. 

99 If there are N numbers of roughness steps, the square of 
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102 where, n d i v = N/L x is the average roughness density 

103 (/cm) at the interface. The other two important scatter- 

104 ing mechanisms are: i) remote charge impurity (RI) scat- 

105 tering, and ii) polar optical phonon scattering. To have 
loe free carriers, the QW is doped remotely. If t is the dis- 
107 tance between the QW and the remotely doped layer, the 
los unscreened matrix element of scattering can be written 
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Polar optical phonon (POP) scattering rates under 
Davydov-Shmushkevich scheme (huo ^> ksT), where 
phonon emission is assumed to be instantaneous, has 
been analytically calculated by Gelmont et. al 14 . The 
only difference between Gelmont' s approach and our cal- 
culation is the form factor, which was calculated by Gel- 
mont et. al using a Fang-Howard type wavefunction, 
whereas, in our case, it is calculated using infinite well 
(hard wall boundary conditions) -type wavefunction. 
Among the scattering potentials considered above, the 
polarization-induced line charge potential and interface- 
roughness potential are anisotropic in nature. This leads 
to anisotropic scattering events which is captured by the 
Kronecker delta function (S q o) appearing in the cor- 
responding matrix elements of scattering. Fig. 2 shows 
a typical anisotropic scattering event where, electron's 
momentum along y direction remains unchanged while 
the momentum along x direction is reversed. For such 
anisotropic scattering, net momentum change is |g| = 
\q x \ = 2k cos#, where 6 is the angle of the incoming wave 
vector with x axis. This is in striking difference with RI 
scattering, where, \q\ = 2/ccos^, where ip — 0' — 0; the 
angle between initial and final wave vector of scattering 
(see inset of Fig. 3b). For line charge scattering and IRs, 
scattering strengths are equally important at all angles 
except at = ±tt/2 where scattering time diverges (no 
scattering), whereas for RI scattering only small angle 
scatterings rates dominate. For such anisotropic scat- 
terers, an angle averaged single relaxation time approx- 
imation (RTA) formalism fails, and one needs to look 
for either variational 16 or numerical solutions of Boltz- 
mann transport equation (BTE). Recently, Schliemann 
and Loss 17 (SL) have proposed an exact solution of BTE 
in the presence of anisotropic scatterers which we use in 
this work. 

To investigate the effect of anisotropic scattering events 
on experimentally measurable transport quantities (such 
as electron mobility, conductivity etc), it is sufficient to 
consider carrier transport along the two principal direc- 
tions x and y. Transport coefficients in any arbitrary di- 
rection, in principle, can be obtained by coordinate trans- 
formation. We first consider charge transport along x. 
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FIG. 3. a): angle (0) dependent scattering time for line-charge 
(black) and IR (blue) scattering, b): RI scattering time as a 
function of angle (ip) between initial and final wave vector of 
scattering. 



153 Under the application of a small electric field E = (£", 0) 

154 in the QW-plane, the non-equilibrium part of carrier dis- 

155 tribution function under the SL scheme can be written 

156 as 



9l(0) 



t)4(*)v( k )-E, 



(6) 



157 where fo is the equilibrium Fermi-Dirac distribution 

158 function, v(k) is the group velocity, and is the kinetic 

159 energy of an electron in the GaN QW. The coefficient 

160 A\(0) is defined as 



i(k) 
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lei where t^t^) is the scattering time parallel (perpendicu- 

162 lar) to applied field calculated using Fermi's golden rule. 

163 Defining the current density as J = 2e/(L x L y ) ^ k gkVk, 

164 the expression for mobility for a degenerate electron gas 

165 along the x direction can be calculated as 
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166 where m* x is the electron effective mass in GaN along 

167 the c axis. A similar expression can be derived for 

168 fiyy. It can be seen that for isotropic scattering, = 

169 A^ = r(kf); which leads to conventional isotropic mo- 

170 bility fi xx = jiyy = er/m*. In general, the integral 

171 appearing in Eq. (9) is evaluated numerically for the 

172 complicated angular dependence of the coefficient A^ (0) . 

173 Nevertheless, analytical expressions of mobilities can be 

174 evaluated under certain approximation. For a very thin 

175 QW, F mn (qa) = 1, and electron mobility for anisotropic 

176 scatterer can be written as 
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FIG. 4. a): IRF scattering vs in-plane correlation length 
A for various carrier densities (Fermi energies), b) electron 
mobility along two principal directions x (blue) and y (black) 
as a function of temperatures. The dashed blue line is the 
longitudinal mobility using variational technique. 
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Where, <?tf is the 2D Thomas- Fermi wave vector 13 , n s 

78 is the equilibrium carrier density in the QW, kp is the 

79 Fermi wave vector, and is the effective Bohr radius 13 . 

The dimensionless integrals ^l^ff^) anc ^ -^(Hf^) can 
be evaluated exactly 18 . 

82 For numerical calculations, a nominal set of parameters 

83 are used to describe the GaN QW roughness: (A, d) = 
(3.19, 5. 2)A, whereas other parameters such as carrier 

85 density (n s ), temperature (T) and the QW thickness (a) 

86 are varied within an experimentally relevant range. An 

87 isotropic effective mass m* x = m* yy = 0.23rao (mo is the 

88 mass of a bare electron) is assumed, since the effective 

89 mass along the c axis {rn\ = 0.228mo) and perpendicular 

to the c axis (m* ± = 0.237mn) has negligible anisotropy 19 
. An important parameter describing the magnitude of 
roughness is the in-plane correlation length A. In Fig. 4, 
we plot the IR scattering rates with in-plane correlation 
length A for different Fermi energies. The scattering rate 
95 exhibits a peak for A « l/k x , where IR scattering ma- 
trix element is maximized, and then decays slowly with 
A on either sides of the peak. For a particular scattering 
rate, it is possible to find two values of A, on each side of 
99 the peak. Since A is a two-valued function of scattering 

200 rates (mobility), it is necessary 20 to determine A from 

201 temperature dependent mobility data. Due to unavail- 

202 ability of transport data for non-polar GaN at present, 

203 we assume A=1.5 nm - a typical value used 15 for polar 

204 GaN HEMTs. It should be noted that the line-charge 

205 scattering rates will have similar but weaker dependence 

206 on A due to inclusion of correlated dipole scattering in 

207 Eq. 4. 

208 The fact that anisotropic scatterers do not hinder elec- 
210 trons moving along y direction, in principle, should result 
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211 in a higher mobility along the y direction. Fig. 4b) cap- 

212 tures this anisotropy. The temperature dependent elec- 

213 tron mobility is evaluated along two principal directions 

214 x and y. At low temperatures, ji yyi limited by RI scatter- 

215 ing, is significantly higher than longitudinal mobility \i xx 

216 which is limited by anisotropic IR and line charge scatter- 

217 ing in addition to RI scattering. At room temperature, 

218 isotropic polar optical phonon scattering tends to reduce 

219 the anisotropy in mobility by equally affecting fi xx and 

220 fiyy. We have assumed distance of remote doping layer 

221 t= 3 nm. Larger values of t will increase the mobility 

222 anisotropy even more at low temperatures by increasing 

223 fiyy exponentially. For example, a value of t = 10 nm 

224 would result in fi yy « 10 5 cm 2 /V.s 21 for n s = 10 12 /cm 2 

225 at low temperatures; approximately 3 times higher mo- 

226 bility compared to t= 3 nm case (see Fig. 4b). The dashed 

227 blue curve in Fig. 3b) shows the mobility (fi xx ) calcu- 

228 lated using the variational principle 16,22 . The difference 

229 of numerical values of mobility between variational and 

230 SL scheme stems out from the fact that variational tech- 

231 niques gives the lower bound of mobility. 

232 Both line-charge and IR scattering matrix elements are 

233 strong decreasing functions of electron's kinetic energy, 

234 which implicitly depends on the magnitude of the mo- 

235 mentum transfer q x in scattering processes. For a degen- 

236 erate electron concentration the mobility is effectively de- 

237 termined the carriers at Fermi level. Hence q x « 2k f = 

238 yj2nn si where n s is the equilibrium carrier density. As a 

239 result, both IR and line charge scattering rates decrease, 

240 which in effect, reduces mobility anisotropy with increas- 

241 ing carrier density n s as illustrated in Fig. 5a). Similarly, 

242 for wide QWs, IR (fim ~ a 6 ) and line charge scattering 

243 are unimportant and polar optical phonon scattering is 

244 the dominant scattering mechanism. Consequently, \i xx 

245 and jiyy tend to approach to the same limiting value 

246 (determined by the POP scattering rate), and mobil- 

247 ity anisotropy is completely washed out (see Fig. 5b)) for 

248 a > 8 nm. At this point, we want to stress upon the fact 

249 that in our numerical calculations, we have used a min- 

250 imal set of parameters for interface roughness (A, d ~ 

251 monolayer). In practice, values of these parameters dif- 

252 fer from sample to sample and experimentally extracted 

253 set of parameters should be used for a more accurate 
25§ quantitative description of transport anisotropy. Two 

256 effects have not been taken account in our model - i) 

257 scattering from charged BSFs and ii) anisotropic strain 

258 at GaN/AIN interface. While scattering from charged 

259 BSFs is expected to enhance transport anisotropy in ad- 

260 dition to the anisotropy presented here, strain-induced 

261 piezoelectric polarization will alter transport anisotropy 

262 by altering bound line charge density at each roughness 

263 step. How to incorporate these two effects in our model 

264 remains an open problem and should be addressed in 

265 future for a more complete and accurate description of 

266 charge transport in non-polar GaN QWs. 

267 In summary, a theory of charge transport in non-polar 

268 GaN QWs has been presented. We show that extended 

269 defects, together with the in-plane polarization of non- 
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FIG. 5. electron mobility along two principal axes x (blue) 
and y (black) — a): as a function of carrier densities and b): 
as a function of well width. Note that mobility anisotropy de- 
creases with increasing carrier density or well width. Dashed 
blue curves represent longitudinal mobility calculated using 
variational technique. 



270 polar GaN-based devices act as anisotropic scattering 

271 centers. At low temperatures, the mobility shows highly 

272 anisotropic behavior for thin QWs. At room temper- 

273 ature, the magnitude of transport anisotropy is reduced 

274 by strong isotropic polar optical phonons scattering. It is 

275 shown that variational technique overestimates the trans- 

276 port anisotropy compared to exact solution of BTE. 
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